Draft version February 1, 2008 

Preprint typeset using LAT^X style emulatcapj v. 04/03/99 



MILLIMETRIC GROUND BASED OBSERVATION OF 
CMBR ANISOTROPY AT <5 = +28° 

G. Romeo 1 , S. Ali 1 , B. Femeni'a, 2 M. Limon, 3 
L. Piccirillo 1 ' 4,5 , R. Rebolo 2 , R. Schaefer 6 
Draft version February 1, 2008 



o 
o 
o 

(N 
> 

o 



> 



O . 
o ■ 

Or 

6 ■ 

a: 



ABSTRACT 

Results from the third campaign of a ground-based multi-band observation of the millimeter emission 
of the sky from Tenerife (Canary Islands) are presented. The instrument consists of a 0.45 meter diameter 
off-axis telescope equipped with a 4-band multi-mode 3 He cooled photometer working at 1.1, 1.3, 2.1 
and 3.1 mm wavelengths. The beam is well approximated by a Gaussian with 1°.35 Full Width Half 
Maximum (FWHM) at all wavelengths. The wide wavelength coverage of our instrument allows us to 
characterize and reduce both the atmospheric and galactic contamination in our data. The CMBR data 
is analyzed in 6 multipole bands whose centers span the range £ = 39 to £ = 134 at the two longest 
wavelengths (2.1 and 3.1 mm). A likelihood analysis indicates that we have detected fluctuations in all 
bands at the two wavelengths. We have evidence of a rise in the angular power spectrum from low £ 
to high I. Our measured spectrum is consistent with current popular theories of large scale structure 
formation, COBE, and other recent balloon-borne experiments with similar wavelength coverage. 

Subject headings: cosmic microwave background — cosmology: observations 



1. INTRODUCTION 



Precise measurements of the angular power spectrum 



Maximum (FWHM) . For details on the instrument se e ( [Fe-| 
mem'a (1998)1) ( |Nicholas (1996)D , ( [Ali et al. (2000)D , and 



of the Cosmic Microwave Background Radiation (CMBR) 
provide strong constraints on cosmological parameters. 
After the era of the d iscov ery of the anisotropy by COBE 
(Smoot et al. (1992) and Bennett ct al. (1996)) we are 
now in the age of spatial spectroscopy of the primordial 
fluctuations. Many experiments have reported detections 
of anisotropy at different angular scales with evidence of 
excess of power at angular scales of about £ ~ 200. 

We designed an experiment sensitive to a broad angular 
range [£ = 39 to 134 ). Our lowest £-band 39±^ par- 
tially overlaps with COBE's highest ^-band. The angular 
range also overlaps with COBE at our longest wavelength 
(3.1 mm). The partial overlap of the parameters of our 
experiment and COBE allows us to perform independent 
consistency checks. 

In this letter, we present the data from our third observ- 
ing campaign of mm- wave CMBR anisotropy from Tenerife 
(Canary Islands). The first campaign in 1993 was intended 
to test the instrumental set-up as well as characterize the 
site. The second campaign in 1994 (see piccirillo ct al 



(1997) and Femema et al. (1998)) measured the anisotropy 



at two £ bands. In the third campaign we extended the 
angular scale coverage from 2 to 6 bands. Observations 
were carried out from May to July, 1996. 

2. INSTRUMENT 

The 1996 Tenerife measurements were performed with 
the same basic telescope of the 1994 campaign. However, 
the beam size has been reduced to 1°.35 Full Width Half 



{Romeo (2000)). Here we briefly summarize its main char- 
acteristics. The sky radiation is collected by the primary 
mirror (0.45 meter diameter). The primary is sinusoidally 
chopped at 2 Hz with a peak-to-peak amplitude of 5°. 7 in 
the sky. The collected radiation is then focused into the 
photometer by means of a fixed off-axis secondary mir- 
ror (0.28 meter diameter). The photometer consists of 
a 4-band He-3 bolometric system operating at 1.1, 1.3, 
2.1 and 3.1 mm wavelengths. The voltage signal from 
each bolometer is sampled synchronously with the sinu- 
soidal motion of the primary mirror. The sampling rate 
is 128 samples/channel/second corresponding to 64 sam- 
ples/channel per sinusoidal cycle of the mirror. In order to 
minimize systematic effects we surrounded the optics with 
two levels of radiation shields. 

3. CALIBRATION 

The primary calibration of the instrument was per- 
formed by observing the mm-wave emission of the Moon. 
Several raster scans of the Moon have been used together 
with a computer simulation to fit the following parame- 
ters: beam size (1°.35), beam throw (5°. 7), the calibra- 
tion constants in Volts/Kelvin for each instrumental band 
and each demodulation, the azimuth and elevation track- 
ing and pointing accuracy. A m odel for the Mo on emission 
temperatures was provided by ( pulkis (1998) ). A correc- 
tion for atmospheric attenuation in our bands was esti- 
mated from skydips as well as from atmospheric models. 
The dominant calibration uncertainty comes from the er- 
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ror in the Moon temperatures (±5%) and the atmospheric 
attenuation (±5%). We add these two errors linearly to 
yield a total calibration uncertainty of ±10%. 

4. OBSERVATIONS 

Two different observing strategies have been used: the 
"Dec40" and the "zenith" modes. The Dec40 mode con- 
sists of observations in drift scan at a fixed elevation cor- 
respo nding to 5 = 40°. This declination has been cxtcn- 



sively | studied at lower frequencies by the Tcncrifc collab 



Sk and the position of the beam in the sky. Observations 
of the Moon signal determine these phases with high accu- 
racy. Each night of observation produces 6 demodulated 
files per bolometric channel. 

The first step in our analysis is to subtract the atmo- 
spheric noise from channels 1, 2, and 3 using the infor- 
mation in channel 4, which is the most sensitive to at- 
mospheric noise. Our technique is essentially the same as 
used in the analysis of our 19 94 c ampaign data (Nicholas 
(1996)|, [Femenia et al. (1998)|, and [Piccirillo et al. (1997)D 



oration (see Gutierrez et al (1999)) as well as from our 
instrument in 1994. For the results of our 1996 campaign 
at declination 40° sec the companion paper ( All et al] 



which w as already an improved version of the technique 
used by Andreani et al. (1991], The det ails o f the pro- 



cedure are to be found in 



Romeo (2000) 



and 



Ali et al 



2000fll~ The data discussed in this Letter are collected (2000)1. A more general theo retical discussion can be found 



during the "zenith" modc:thc telescope is fixed to observe 
the zenith and the sinusoidal chopping throws the beam 
along the North-South direction. The East- West sky ro- 
tation provides the sampling of a strip of the sky defined 
by Oh < RA < 2Ah. This obse rving strategy is similar t o 
the Saskatoon experiment (see Netterfield et al. (1997) ), 
although we observe at a different declination. 

Data were collected from May through June 1996 when 
atmospheric fluctuations are at their minimum. Any fur- 
ther contamination was minimized by using only the data 
collected at night, when the Sun was well below the hori- 
zon. About 241 hours of data (from 23 nights of observa- 
tion), corresponding to about 1000 square degrees of sky 
observed at 4 different wavelengths were collected by using 
the zenith- mode. 

Ou r zenith observing strategy has the advantage that 



the t qlescope does not move during the CMBR observa- 
tions. This virtually eliminates any changes in the sidelobc 
pick-up that can contaminate the data due to the motion 
of the telescope. 

5. DATA REDUCTION AND ANALYSIS 

Data are demodulated using the following set of orthog- 
onal functions 

L n (ujt) = cos{nujt) + i sm(nut) = e {inut) (1) 

where u> is the chopping frequency. The output of the 
n — th synthesized beam is: 



fc=l 



(2) 



where Sk are the discrete samples of the bolometer volt- 
ages, N=64 is the number of samples per sinusoidal cy- 
cle and n=1...6 are the 6 demodulation vectors indices, 
which can be used to reconstruct a 2-dimensional map of 
anisotropy ( Romeo et al. (2000) ) . Our choice of orthogo- 
nal functions allows us to express the demodulated output 
according to: 
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where the ipn = tan 1 (b n /a n ) are phases which account 
for the instrumental delays between the time ordered data 



in |Mclchiorri et al. (1996) . The atmospheric subtraction 
can only be applied when the atmospheric conditions are 
stable and there is a strong correlation (> 75%) between 
channel 4 and the other channels. In the lower harmonics, 
the power spectra of the atmospheric subtracted data show 
a strong improvement in the signal to noise ratio. For the 
higher harmonics (n = 4, 5, and 6), we did not need the 
atmospheric subtraction because the high degree of spa- 
tial differencing breaks up the atmospheric noise, whose 
power decreases sharply with scale. (Atmospheric noise 
follows a Kolmogorov power spectrum which decreases as 
the wavenumber to the -8/3 power. ) A more in-depth 
analysis of the a tmospheric nois e will be discusse d in even 



greater detail in Ali (2000) and Ali et al. (2000) 



Long term drifts in the data are then remove d by fit- 
ting a combin atio n of sine and co sine functions (Femenia 



et al. (1998) and Femenia (1998)). The spatial frequen- 



cies used correspond to angles of 90, 60, 30, 30, 30 30 
degrees R. A. for demodulations 1-6 respectively, angular 
scales where the window functions are negligible. Data 
are finally stacked together and binned in RA intervals of 
about 1°.4. 

The final data set is analyzed by performing a likeli- 
hood analysis of each individual demodulation for the two 
longest wavelength bands (3.1 and 2.1 mm) corresponding 
to 6 £ ranges centered from £=39 to ^=134. In fig. l] we 
show the first demodulation of the bands together with the 
demodulations obtained by simulating our observing strat- 
egy on the DIRBE 240 \im map . The feature visible in 
the DIRBE demodulations at RA about 20h is the region 
corresponding to the crossing of the galactic plane. We 
see that the same feature is visible in the 1.3 mm band 
and somewhat less visible with increasing wavelength to 
2.1 and 3.1 mm. 

We can use the ratios of intensities I (DIRBE) /I (mm) to 
extrapolate the galactic signal in the region of 12 < RA < 
19 where we analyze the data for CMBR anisotropy. The 
result of the extrapolation is shown in fig. ||. We see that, 
at high galactic latitude, in the two longest wavelength 
bands (1 and 2) the galactic contamination should be neg- 
ligible while the band at 1.3 mm might be contaminated 
by residual galactic dust emission and is not included in 
this analysis. 

Computing the various window functions is the last step 
needed to perform the likelihood analysis. The relative 
simplicity in the observing strategy comes at the cost of 
some algebra in computing the window functions for each 
individual demodulation at any angular lag. Following 
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( [White and Srcdnicki (1994) ), the expression for the win- 
dow function of each demodulation is: 



2 r>2 ) n» 



(£ + m)\ 



cos(m(0 4 —</>j)) 



sin(^) 

mA0 



(4) 



where n is the integer for the corresponding demodulation, 
ki and kj are respectively the unit vectors identifying right 
ascension bin i and j, N is the norm alization determined 
in the same way as for Saskatoon [see White and Srednickj 
(1994F 



exp 



-e(£+l)a 2 



is the beam profile function with <t—0.°57, our Gaussian 
beam width, and 



Q 7 £.m(@o) — 



2tt 
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L n (t)PP[cos(6(t))]dt m = 0...£ 



wher e L n (t) is the lock-in function in equation [T| 



P/ n [cos(6*(i))] is the associated Legendre polynomials, A<f> 
is the binning size in Right Ascension, and finally 6(t) = 
9 + asia.(wt), with O = 28°, and a = 2°.84. 

The window functions are estimated numerically. 

6. STATISTICAL ANALYSIS 

The statistical analysis used to determine the amplitude 
of the fluctuation in the CMB was performed by using the 
Likelihood defined as follows: 
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where t is the vector containing the demodulated and 
stacked data. We maximize the likelihood L for channels 
1 and 2 both individually and collectively. Each channel 
has 72 1.4 degree bins of stacked data for each harmonic, 
so the covariance matrix M is a 72 x 72 (144 x 144) matrix 
for the individual (combined) likelihood, respectively. The 
combined channel analysis includes terms for the correla- 
tion between channels. The covariance matrix is the sum 
of two matrices, signal ($y) and noise (Nij). The signal 
is given by 



where WP(fti ■ n.2) is the window function, and Cg are the 
usual 



the signal autocorrelation function, we used a flat spa- 
tial temperature fluctuation spectrum parametrized by the 
quadru pole amplitude AT (see e. g., schaefer and deLaix 



(1996)| ) 



6x/W (AT)2 21 
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AT is then fit to the data and the result is converted to 
band power. We obtain very consistent results separately 
for channels 1 and 2 which are shown in table 1 along 
with the combined channel results. Results for channel 1 
and 2 are in Table The reported 68% co nfidence levels 



have b een calculated in the usual way (see Church et al 



1997)), and do not include calibration uncertainty. In Fig- 



ure |3j, we show two theoretical structure formation mod- 
els for comparison, the "standard" cold dark model and a 
spatially flat cosmological constant (A) dominated model. 
We see that our data are consistent with the rise to the 
"Doppler peak" expected in currently popular adiabatic 
theories of structure formation (like the models shown.) 

The cosmic origin of our signal has been further checked 
with the following null test: we divided the set of the 
stacked data files in two subsets; the likelihood analysis 
of their sum and difference produced respectively results 
consist ent with those in Table [j], and with zero (see Romcc 



(2000)) . 



We can also calculate the overlap window functions be- 
tween different lock-in patterns and angular lags and use it 
to estimate the correlations between harmonics, similar to 



tumpuiatmc augulai 



jflicieiiLs. — Tfo 



the method in ( Ncttcrficld (1995) ). The amplitude of the 
correlated harmonics was calculated using these overlap 
window functions and a flat CMB spectrum, normalized 
so the diagonal elements are unity. We found that correla- 
tions between signals measured by the various harmonics 
are for the most part negligible (< 1%). Even-odd pairs 
of harmonics are completely uncorrelated. There is some 
non-negligible correlation between pairs of either odd-odd 
or even-even harmonics. The correlations between har- 
monic pairs 1-3, 2-4, 3-5, and 4-6 are 0.42, 0.53, 0.61 and 
0.67, respectively. We verified that these correlation terms 
are detectable in a likelihood analysis of the data and are 
consistent with our estimates. The correlations between 
pairs 1-5 and 2-6 are much smaller: 0.03 and 0.04, re- 
spectively. If one wants to be conservative in using these 
points to constrain cosmological parameters, one can use 
one even-odd pair of harmonics and be assured that those 
two points will be completely independent. 

Finally we compare our results with other experiments 
in Figure 3. We see that our lowest 1-band temperature 
measurement is completely consistent with COBEs high- 
est 1-band temperature fluctuation measurement. Fur- 
thermore, we plot results from two recent balloon ex- 
pcrimcnts done at similar w avelengths, BOO MERanG 
flde Bernardis, et al. (2000)| ), and MAXIMA ( [Hanany 



noise covariance matrix is 



et al. (2000) J. (For a more complete summary of ex- 



per imental results, see Max Tegmark's CMB web sit e 



(TiUjCij 



where is the correlation function obtained from the 
data. Because of the symmetry of the covariance ma- 
trix we can use Cholesky factorization to invert M. For 



at http:/ /www. sns.ias.edu/max/cmb/experiments. html.) 
We see that our results are remarkably consistent with 
those of other experiments, which observed a different part 
of the sky, and used completely different experimental and 
analysis techniques. 
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7. CONCLUSION 

Improvements in our instrumental setup over our 1994 
campaign, as well as refinements in our data analysis pro- 
cedure, have allowed us to detect temperature fluctuations 
over a wide range of angular scales. We have used data 
from our 4-band detector to reduce the atmospheric noise 
contamination and to delineate regions of our scans for 
which the galactic contamination is insignificant. We find 
positive detections of residual temperature fluctuations in 
6 angular ranges in our two longest wavelength bands (2.1 
mm and 3.1 mm). These measurements, summarized in 



table 1, are our main result. The derived temperatures 
are consistent with those seen in other experiments . Our 
results essentially show a monotonically increasing ampli- 
tude of temperature fluctuation with 1, and are consistent 
with some currently popular adiabatic theories of cosmo- 
logical structure formation. 
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staff of the Observatorio del Teide, and finally R. Hoyland, 
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Fig. 1. — The first harmonics of thermodynamic temperature differences in the three channels. The thick line is the result of simulating our 
observing pattern on the DIRBE 240 /i m map of the sky, scaled to the corresponding wavelength. The scaling factors are respectively 5.0, 
5.4, and 30.8 /iK/MJy/sr for our 3.1 mm, 2.1mm, and 1.3 mm wavelength channels. 

Fig. 2. — Comparison between Galactic Dust Emission and CMB anisotropy. The galactic dust emission for each channel has been evaluated 
in the galactic plane crossing region (RA = 19-21 hours) and extrapolated to the other region (RA = 12-19 hours), using information from 
the first harmonic on real data and the simulated observations on the DIRBE map (see figure 1). This is done for channels 1, 2, and 3. The 
spectral index of the galactic dust emission is estimated to be 1.5 ± 1. from fitting a power law line to the data shown above. 



FlG. 3. — Band power vs. 1 for all harmonics. We have obtained the theoretical curves using CMBFAST (3eljak and Zaldarriaga (1996)). 
Shown are "standard CDM" (n CDM , tt A , f2 6 , n, h = Hn/-\() rik m .9- 1 Mpr' 1 ) = (0 95, 0.0 (1 0.05, 1.0, 0.65) and A CDM (0.35, 0.60, 0.05, 
1.0, 0.65). The data shown are-Emm experiment, r^npF. qTcgmark ancLHaxajltog ( 1QQ7 ) ). ar y* two other recent experiments using similar 
wavelengths - BOOMERanG ( ie Bernardis, ct al. (2000)[) and MAXIMA (Hanany ct al. (2000)). 
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Table 1 

Summary of Band Powers for channel 1 and 2 as determined by the likelihood analysis. 



Harmonic i Channel 1 Channel 2 Channel 1 + Channel 2 
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■ Galactic Emission at 19 < ra < 21 
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